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ABSTRACT 

Use of the highly sensitive Hokupa'a/Gemini curvature wavefront sensor has allowed direct adaptive 
optics (AO) guiding on very low mass (VLM) stars with SpT=M8.0-L0.5. A survey of 39 such objects 
detected 9 VLM binaries (7 of which were discovered for the first time to be binaries). Most of these 
systems are tight (separation < 5 AU) and have similar masses (AKs < 0.8 mag; 0.85 < q < 1.0). 
However, 2 systems (LHS 2397a and 2M2331016-040618) have large AKs > 2.4 mag and consist of a 
VLM star orbited by a much cooler L7-L8 brown dwarf companion. Based on this flux limited (Ks < 12 
mag) survey of 39 M8.0-L0.5 stars (mainly from the 2MASS sample of Gizis et al. (2000)) we find a 
sensitivity corrected binary fraction in the range 15 ±7% for M8.0-L0.5 stars with separations > 2.6 AU. 
This is slightly less than the 32 ± 9% measured for more massive M0-M4 dwarfs over the same separation 
range (Fischer & Marcy 1992). It appears M8.0-L0.5 binaries (as well as L and T dwarf binaries) have 
a much smaller semi-major axis distribution peak (~ 4 AU) compared to more massive M and G dwarfs 
which have a broad peak at larger ~ 30 AU separations. We also find no VLM binary systems (defined 
here as systems with M to t < 0.185M Q ) with separations > 15 AU. 

We briefly explore possible reasons why VLM binaries are slightly less common, nearly equal mass, 
and much more tightly bound compared to more massive binaries. We investigate the hypothesis that 
the lack of wide (a > 20 AU) VLM/brown dwarf binaries may be explained if the binary components 
were given a significant differential velocity kick. Such a velocity kick is predicted by current "ejection" 
theories, where brown dwarfs are formed because they are ejected from their embryonic mini-cluster and 
therefore starved of accretion material. We find that a kick from a close triple or quadruple encounter 
(imparting a differential kick of ^ 3 km/s between the members of an escaping binary) could reproduce 
the observed cut-off in the semi-major axis distribution at ~ 20 AU. However, the estimated binarity 
(< 5%; Bate et al. (2002)) produced by such ejection scenarios is below the 15 ± 7% observed. Similarly, 
VLM binaries could be the final hardened binaries produced when a mini-cluster decays. However, the 
models of Sterzik & Durisen (1998); Durisen, Sterzik, & Pickett (2001) also cannot produce a VLM 
binary fraction above ~ 5%. The observed VLM binary frequency could possibly be produced by cloud 
core fragmentation. Although, our estimate of a fragmentation-produced VLM binary semi-major axis 
distribution contains a significant fraction of "wide" VLM binaries with a > 20 AU in contrast to 
observation. In summary, more detailed theoretical work will be needed to explain these interesting 
results which show VLM binaries to be a significantly different population from more massive M & G 
dwarf binaries. 

Subject headings: instrumentation: adaptive optics — binaries: general — stars: evolution — stars: 
formation — stars: low-mass, brown dwarfs 



1. INTRODUCTION 

Since the discovery of Gl 229B by Nakajima et al. (1995) 
there has been intense interest in the direct detection of 
brown dwarfs and very low mass (VLM) stars and their 
companions. According to the current models of Burrows 
et al. (2000) and Chabrier et al. (2000), stars with spectral 
types of M8.0-L0.5 will be just above the stellar /substellar 
boundary. However, modestly fainter companions to such 
primaries could themselves be substellar. Therefore, a sur- 
vey of M8.0-L0.5 stars should detect binary systems con- 
sisting of VLM primaries with VLM or brown dwarf sec- 
ondaries. 

The binary frequency of M8.0-L0.5 stars is interesting 
in its own right since little is known about how common 
M8.0-L0.5 binary systems are. It is not clear currently if 
the M8.0-L0.5 binary separation distribution is similar to 
that of M0-M4 stars; in fact, there is emerging evidence 
that very low mass L & T dwarf binaries tend to have 



smaller separations and possibly lower binary frequencies 
compared to more massive M and G stars (Martin, Brand- 
ner, & Basri 1999; Reid et al. 2001a; Burgasser et al. 2003). 

Despite the strong interest in such very low mass (VLM) 
binaries (Mtot < O.185Af ), only 24 such systems are 
known (see Table 4 for a complete list). A brief overview of 
these systems starts with the first double L dwarf system 
which was imaged by HST/NICMOS by Martin, Brand- 
ner, & Basri (1999). A young spectroscopic binary brown 
dwarf (PPL 15) was detected in the Pleiades (Basri & 
Martin (1999)) but this spectroscopic system is too tight 
to get separate luminosities for each component. A large 
HST/NICMOS imaging survey by Martin et al. (2000) of 
VLM dwarfs in the Pleiades failed to detect any brown 
dwarf binaries with separations > 0.2" (> 27 AU). Detec- 
tions of nearby field binary systems were more successful. 
The nearby object Gl 569B was resolved into a 0.1" (1 AU) 
binary brown dwarf at Keck and the 6.5m MMT (Martin 
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et al. (1999); Kenworthy et al. (2001); Lane et al. (2001)). 
Keck seeing-limited NIR imaging marginally resolved two 
more binary L stars (Koerner et al. (1999)). A survey with 
WFPC2 detected four more (three newly discovered and 
one confirmed from Koerner et al. (1999)) tight equal mag- 
nitude binaries out of a sample of 20 L dwarfs (Reid et al. 
(2001a)). From the same survey Reid et al. (2002) found 
a M8 binary (2M1047; later discovered independently by 
our survey). Guiding on HD 130948 with adaptive optics 
(AO), Potter et al. (2002a) discovered a companion binary 
brown dwarf system. Recently, Burgasser et al. (2003) 
have detected two T dwarf binaries with HST. Finally, 12 
more L dwarf binaries have been found by analyzing all 
the currently remaining HST/WFPC2 data collected on 
L dwarfs (Bouy et al. 2003). Hence, the total number of 
binary VLM stars and brown dwarfs currently known is 
just 24. Of these, all but one have luminosities known for 
each component, since one is a spectroscopic binary. 

Here we have carried out our own high-spatial resolu- 
tion binary survey of VLM stars employing adaptive op- 
tics. In total, we have detected 12 VLM binaries (10 of 
these are new discoveries) in our AO survey of 69 VLM 
stars. Three of these systems (LP415-20, LP475-855, 
& 2MASSW J1750129+442404) have M7.0-M7.5 spectral 
types and are discussed in detail elsewhere (Siegler et al. 
2003). In this paper, we discuss the remaining 9 cooler 
binaries with M8.0-L0.5 primaries detected in our survey 
(referred herein as M8.0-L0.5 binaries even though they 
may contain L1-L7.5 companions; see Table 2 for a com- 
plete list of these systems). 

Two of these systems (2MASSW J0746425+200032 
and 2MASSW J1047127+402644) were in our sample 
but were previously imaged in the visible by HST 
and found to be binaries (Reid et al. (2001a, 2002)). 
Here we present the first resolved IR observations of 
these two systems and new astrometry. The seven 
remaining systems were all discovered to be bina- 
ries during this survey. The first four systems dis- 
covered in our survey (2MASSW J1426316+155701, 
2MASSW J2140293+162518, 2MASSW J2206228-204705, 
and 2MASSW J2331016-040618) have brief descriptions 
in Close et al. (2002b). However, we have re-analyzed 
the data from Close et al. (2002b) and include it here 
for completeness with slightly revised mass estimates. 
The very interesting M8/L7.5 system LHS 2397a discov- 
ered during this survey is discussed in detail elsewhere 
(Freed, Close, & Siegler 2003) yet is included here for 
completeness. The newly discovered binaries 2MASSW 
J1127534+741107 and 2MASSW J1311391+803222 are 
presented here for the first time. 

These nine M8.0-L0.5 binaries are a significant addi- 
tion to the other very low mass M8-T6 binaries known 
to date listed in Table 4 (Basri & Martin 1999; Martin, 
Brandncr, & Basri 1999; Koerner et al. 1999; Reid et al. 
2001a; Lane et al. 2001; Potter et al. 2002a; Burgasser et 
al. 2003; Bouy et al. 2003). With relatively short periods 
our new systems will likely play a significant role in the 
mass-age-luminosity calibration for VLM stars and brown 
dwarfs. It is also noteworthy that we can start to charac- 
terize this new population of M8.0-L0.5 binaries. We will 
outline how VLM binaries are different from their more 
massive M and G counterparts. Since VLM binaries are 



so tightly bound we hypothesize that very little dynami- 
cal evolution of the distribution has occurred since their 
formation. We attempt to constrain the formation mech- 
anism of VLM stars and brown dwarfs from our observed 
semi-major axis distribution and binarity statistics. 

2. AN AO SURVEY OF NEARBY M8.0-L0.5 FIELD STARS 

As outlined in detail in Close et al. (2002a), we utilized 
the University of Hawaii curvature adaptive optics system 
Hokupa'a (Graves et al. 1998; Close et al. 1998), which was 
a visitor AO instrument on the Gemini North Telescope. 
This highly sensitive curvature AO system is well suited 
to locking onto nearby, faint (V ~ 20), red (V — I > 3), 
M8.0-L0.5 stars to produce ~ 0.1" images (which are close 
to the 0.07" diffraction- limit in the K' band). We can 
guide on such faint (/ ~ 17) targets with a curvature AO 
system (such as Hokupa'a) by utilizing its zero-read noise 
wavefront sensor (for a detailed explanation of how this is 
possible see Siegler, Close, & Freed (2002)). We utilized 
this unique capability to survey the nearest extreme M 
and L stars (M8.0-L0.5) to characterize the nearby VLM 
binary population. 

Here we report the results of all our Gemini observing 
runs in 2001 and 2002. We have observed all 32 M8.0- 
M9.5 stars with Ks < 12 mag from the list of Gizis et al. 
(2000). It should be noted that the M8.0-M9.5 list of Gizis 
et al. (2000) has some selection constraints: galactic lati- 
tudes are all > 20 degrees; and from < RA < 4.5 hours 
DEC < 30 degrees; and there are gaps in the coverage 
due to the past availability of the 2MASS scans. A bright 
L0.5 dwarf with Ks < 12 was also observed (selected from 
Kirkpatrick et al. (2000)). Six additional bright (Ks < 12) 
M8.0-M9.5 stars were selected from Reid et al. (2002) and 
Cruz et al. (2003). In total 39 M8.0-L0.5 stars have now 
been imaged at high resolution (~ 0.1") with AO compen- 
sation in our survey For a complete list of these M8-L0.5 
target stars see Table 1 (single stars) and Table 2 (stars 
found to be binaries). 

Nine of our 39 targets were clearly tight binaries (sep < 
0.5"). We observed each of these objects by dithering 
over 4 different positions on the QUIRC 1024 x 1024 NIR 
(1 — 2.5jum) detector (Hodapp et al. 1996) which has a 
0.0199"/pixcl plate scale at Gemini North. At each po- 
sition we took 3x10s exposures at J, H, K' , and 3x60s 
exposures at H, resulting in unsaturated 120s exposures 
at J, H, and K' with a deep 720s exposure at H band for 
each binary system. 

3. REDUCTIONS 

We have developed an AO data reduction pipeline in the 
IRAF language which maximizes sensitivity and image res- 
olution. This pipeline is standard IR AO data reduction 
and is described in detail in Close et al. (2002a). 

Unlike conventional NIR observing at an alt-az telescope 
(like Gemini), we disable the Cassegrain rotator so that 
the pupil image is fixed w.r.t. the detector. Hence the 
optics are not rotating w.r.t. the camera or detector. In 
this manner the residual static PSF aberrations are fixed 
in the AO images enabling quick identification of real com- 
panions as compared to PSF artifacts. The pipeline cross- 
correlates and aligns each image, then rotates each image 
so north is up and east is to the left, then median combines 
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the data with an average sigma clip rejection at the ±2.5c 
level. By use of a cubic-spline interpolator the script pre- 
serves image resolution to the < 0.02 pixel level. Next 
the custom IRAF script produces two final output images, 
one that combines all the images taken and another where 
only the sharpest 50% of the images are combined. 

This pipeline produced final unsaturated 120s expo- 
sures at J (FWHM ~ 0.15") , H (FWHM ~ 0.14"), 
and K' (FWHM ~ 0.13") with a deep 720s exposure 
(FWHM ~ 0.14") at H band for each binary system. The 
dithering produces a final image of 30x30" with the most 
sensitive region (10 x 10") centered on the binary. Figures 
1 and 2 illustrates K' images of each of the systems. 

In Table 2 we present the analysis of the images taken 
of the 9 new binaries from our Gemini observing runs. 
The photometry was based on DAOPHOT PSF fitting 
photometry (Stetson 1987). The PSFs used were the re- 
duced 12x10s unsaturated data from the next (and pre- 
vious) single VLM stars observed after (and before) each 
binary. The PSF stars always had a similar IR brightness, 
a late M spectral type, and were observed at a similar 
airmass. The resulting Amagnitudes between the compo- 
nents are listed in Table 2; their errors in Amag are the 
differences in the photometry between two similar PSF 
stars. The individual fluxes were calculated from the flux 
ratio measured by DAOPHOT. We made the assumption 
AK' ~ AKs, which is correct to 0.02 mag according the 
models of Chabrier et al. (2000). Assuming AK' = AKs 
allows us to use the 2MASS integrated Ks fluxes of the 
blended binaries to solve for the individual Ks fluxes of 
each component (see Table 3). 

The platescale and orientation of QUIRC was deter- 
mined from a short exposure of the Trapezium cluster in 
Orion and compared to published positions as in Simon, 
Close, & Beck (1999). From these observations a platescale 
of 0.0199 ± 0.0002"/pix and an orientation of the Y-axis 
(0.3 ± 0.3 degrees E of north) was determined. Astrom- 
etry for each binary was based on the PSF fitting. The 
astrometric errors were based on the range of the 3 values 
observed at J, H, and K' and the systematic errors in the 
calibration added in quadrature. 

4. ANALYSIS 

4.1. Are the companions physically related to primaries? 

Since Gizis et al. (2000) only selected objects > 20 de- 
grees above the galactic plane, we do not expect many 
background late M or L stars in our images. In the 3.6xl0 4 
square arcsecs already surveyed, we have not detected a 
very red J — Ks > 0.8 mag background object in any 
of the fields. Therefore, we estimate the probability of a 
chance projection of such a red object within < 0.5" of the 
primary to be < 2xl0~°. Moreover, given the rather low 
space density (0.0057 ± 0.0025pc~ 3 ) of L dwarfs (Gizis et 
al. 2001), the probability of a background L dwarf within 
< 0.5" of any of our targets is < 10~ 16 . We conclude that 
all these very red, cool objects are physically related to 
their primaries. 

4.2. What are the distances to the binaries? 

Unfortunately, there are no published trigonometric par- 
allaxes for six of the nine systems. The three systems with 
parallaxes are: 2M0746 (Dahn et al. 2002); LHS 2397a 



(van Altena, Lee, & Hofileit 1995); and 2M2331 (associ- 
ated with a Hipparcos star HD221356 Gizis et al. (2000)). 
For the remaining six systems without parallaxes we can 
estimate the distance based on the trigonometric paral- 
laxes of other well studied M8.0-L0.5 stars from Dahn et 
al. (2002). The distances of all the primaries were deter- 
mined from the absolute Ks magnitudes (using available 
2MASS photometry for each star with trigonometric par- 
allaxes from Dahn et al. (2002)), which can be estimated 
by M Ks = 7.71 + 2.U( J-Ks) for M8.0-L0.5 stars (Siegler 
et al. 2003). This relationship has a la error of 0.33 mag 
which has been added in quadrature to the J and Ks pho- 
tometric errors to yield the primary component's Mk s val- 
ues in Table 3 and plotted as crosses in Figures 3 - 11. As 
can be seen from Table 3 all but one of our systems are 
within 29 pc (the exception is 2M1127 at ~ 33 pc). 

4.3. What are the spectral types of the components? 

We do not have spatially resolved spectra of both com- 
ponents in any of these systems; consequently we can 
only try to fit the Mk s values in Table 3 to the rela- 
tion SpT = 3.97M Ks - 31.67 which is derived from the 
datasct of Dahn et al. (2002) by Siegler et al. (2003). 
Unfortunately, the exact relationship between Mk.s and 
VLM/brown dwarf spectral types is still under study. It 
is important to note that these spectral types are only a 
guide since the conversion from Mk s to spectral type car- 
ries at least ±1.5 spectral subclasses of uncertainty. For- 
tunately, none of the following analysis is dependent on 
these spectral type estimates. 

It is interesting to note that six of these secondaries are 
likely L dwarfs. In particular 2M2331B is likely a L7 and 
LHS 2397aB is likely a L7.5. Both 2M2331B and LHS 
2397aB are very cool, late L companions. 

4.4. What are ages of the systems? 

Estimating the exact age for any of these systems is dif- 
ficult since there are no Li measurements yet published 
(which could place an upper limit on the ages). An ex- 
ception to this is LHS 2397a for which no Li was detected 
Martin, Rcbolo, & Magazzu (1994). For a detailed discus- 
sion on the age of LHS 2397a see Freed, Close, & Siegler 
(2003). For each of the remaining systems we have conser- 
vatively assumed that the whole range of common ages in 
the solar neighborhood (0.6 to 7.5 Gyr) may apply to each 
system (Caloi et al. 1999). However, Gizis et al. (2000) 
observed very low proper motion (Vt an < 10 km/s) for the 
2M1127, 2M2140, and 2M2206 systems. These three sys- 
tems are among the lowest velocity M8's in the entire sur- 
vey of Gizis et al. (2000) suggesting a somewhat younger 
age since these systems have not yet developed a signifi- 
cant random velocity like the other older (~ 5 Gyr) M8.0- 
L0.5 stars in the survey. Therefore, we assign a slightly 
younger age of 3. 0+^4 Gyr to these 3 systems, but leave 
large error bars allowing ages from 0.6-7.5 Gyr (~ 3 Gyr is 
the maximum age for the kincmatically young stars found 
by Caloi et al. (1999)). The other binary systems 2M0746, 
2M1047, 2M1311, and 2M2331 appear to have normal V tan 
and are more likely to be older systems. Hence we assign 
an age of 5.0jl|;| Gyr to these older systems (Caloi et al. 
1999). It should be noted that there is little significant 
difference between the evolutionary tracks for ages 1 — 10 
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Gyr when SpT < LO (Chabrier et al. 2000). Therefore, 
the exact age is not absolutely critical to estimating the 
approximate masses for M8.0-L0.5 stars (see Figure 3). 

4.5. The masses of the components 

To estimate masses for these objects we will need to 
rely on theoretical evolutionary tracks for VLM stars and 
brown dwarfs. Calibrated theoretical evolutionary tracks 
are required for objects in the temperature range 1400- 
2600 K. Recently such a calibration has been performed 
by two groups using dynamical measurements of the M8.5 
G1569B brown dwarf binary. From the dynamical mass 
measurements of the G1569B binary brown dwarf (Ken- 
worthy et al. 2001; Lane et al. 2001) it was found that 
the Chabrier et al. (2000) and Burrows et al. (2000) evo- 
lutionary models were in reasonably good agreement with 
observation. In Figures 3 to 11 we plot the latest DUSTY 
models from Chabrier et al. (2000) which have been spe- 
cially integrated across the Ks filter so as to allow a direct 
comparison to the 2MASS Ks photometry (this avoids the 
additional error of converting from Ks to K for very red 
objects). We extrapolated the isochrones from 0.10 to 0.11 
Mq to cover the extreme upper limits of some of the pri- 
mary masses in the figures. 

We estimate the masses of the components based on the 
age range of 0.6-7.5 Gyr and the range of Mk s values. The 
maximum mass relates to the minimum Mk s and the max- 
imum age of 7.5 Gyr. The minimum mass relates to the 
maximum Mk s and the minimum age of 0.6 Gyr. These 
masses are listed in Table 3 and illustrated in Figures 3 to 
11 as filled polygons. 

At the younger ages (< lGyr), the primaries may be on 
the stellar/substellar boundary, but they are most likely 
VLM stars. The substellar nature of the companion is very 
likely in the case of 2M2331B and LHS 2397aB, possible in 
the cases of 2M0746B, 2M1426B, and 2M2140B, and un- 
likely in the cases of 2M1047B, 2M1127B, 2M1311B, and 
2M2206B which all appear to be VLM stars like their pri- 
maries. Hence two of the companions are brown dwarfs, 
three others may also be substellar, and four are likely 
VLM stars. 

5. DISCUSSION 

5.1. The binary frequency of M8.0-L0.5 stars 

We have carried out the largest flux limited (Ks < 12) 
high spatial resolution survey of M8.0-L0.5 primaries. 
Around these 39 M8.0-L0.5 targets we have detected 9 
systems that have companions. Since our survey is flux 
limited we need to correct for our bias toward detect- 
ing equal magnitude binaries that "leak" into our sample 
from further distances. For example, an equal magnitude 
M8 binary could have an integrated 2MASS magnitude 
of Ks=12 mag but be actually located at 36 pc whereas 
a single M8 star of Ks=12 would be located just 26 pc 
distant. Hence our selection of Ks < 12 leads to incom- 
pleteness of single stars and low mass ratio (q=M 2 /Mi) 
binaries past D ~ 26 pc. More exactly, 88% of our bi- 
nary systems are within 29.1 pc (distances calculated using 
only the primary's apparent magnitude) and 88% of our 
single stars are within 23.4 pc. Therefore, we are prob- 
ing ~ (29. 1/23. 4) 3 = 1.92 times more volume with the 
brighter (combined light) of the binaries compared to the 



single (hence fainter) M8.0-L0.5 stars. Hence, the cor- 
rected binary frequency is 9/39/1.92 = 12 ± 4% (where 
the error is only Poisson error). 

There is another selection effect due to the instrumen- 
tal PSF which prevents detection of very faint compan- 
ions very close to the primaries. At the smallest sepa- 
rations of 0.1 — 0.2" we are only sensitive to relatively 
bright companions of AK' < lmag. Much fainter com- 
panions (AK' <~ bmag) can be detected at slightly wider 
(~ 0.25") separations, and very low mass companions 
(AH <~ lOmag) could be detected at ~ 1" separations 
in our deep 720s H images. Therefore, we are likely in- 
sensitive to faint (AK' > 1.0) companions in the sepa- 
ration range of 0.1 — 0.2". However, if we assume that 
the mass ratio (q) distribution for M8.0-L0.5 stars is close 
to flat (as it is for M0-M4 binaries; Fischer & Marcy 
(1992)), then we would expect at least as many binaries 
with AK' > 1.0 as AK' < 1.0 mag. Although we do 
not have enough data currently to definitely derive the 
q distribution for M8.0-L0.5 binaries, we can note that 
for the four systems with separations > 0.2" we observed 
an equal number of AK' > 1.0 as AK' < 1.0 mag sys- 
tems. So it appears reasonable that there should also be 
an equal number of AK' > 1.0 as AK' < 1.0 mag sys- 
tems in the range 0.1 — 0.2". Consequently, based on our 
detection of five systems with AK' < 1 mag with sepa- 
rations of 0.1 — 0.2" we would expect to have ~ 5 sys- 
tems with AK' > 1.0 in the range 0.1 — 0.2". In reality 
we detected no systems with AK' > 1.0 with separations 
0.1 — 0.2". Therefore, to correct for instrumental insensi- 
tivity we need to increase the number of binary systems by 
5 in the range 0.1 — 0.2". Based on this assumption about 
the mass ratio distribution there should be ~ 10 binaries 
from 0.1 — 0.2" when correcting for our instrumental in- 
sensitivity. Therefore, the total count for all separations 
> 0.1" should be 14 ± 4 systems assuming a Poisson er- 
ror. Therefore, the corrected M8.0-L0.5 binary frequency 
would be 14/39/1.92 = 19 ± 7% for separations > 0.10" 
or > 2.6 AU. Hence we have a range of possible volume- 
limited binary frequencies (BF) from 12% (q ~ 1) up to 
19% (where q is assumed to be flat necessitating a large 
correction for inscnsitivity). 

As a check we can re-derive these values of the true 
volume-limited BF in a manner after Burgasser et al. 
(2003) for the two limiting cases of q = 1 up to a flat 
q distribution. We can write, 

BF = -BF obs /(BF obs + a(BF obs - 1)) (1) 

where BF obs is the observed uncorrected binarity 
(BF obs = 9/39 = 23%), and where a is the fractional 
increase in volume sampled for binaries with a flux ratio 
p = Fb/Fa and a flux ratio distribution f(p). As in (Bur- 
gasser et al. 2003) we can calculate a by, 

Jo f(P) d P 

We consider two limiting cases for the f(p) distribution: 
1) if all the systems are equal magnitude (q = p = 1) 
then a = 2^/ 2 ) = 2.8 and the BF=12%; 2) if there is 
a flat f(p) distribution then a — 1.9 and the BF=19%. 
Consequently, the binary frequency range is 12-19% which 
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is identical to the range estimated above. Later we will 
see (Figure 14) the true f(p) is indeed a compromise be- 
tween flat and unity; hence we split the difference and 
adopt a binary frequency of 15 ± 7% where the error is 
the Poisson error (5%) added in quadrature to the (<~ 4%) 
uncertainty due to the possible range of the q distribution 
(1.9 < a < 2.8). It appears that for systems with separa- 
tions 2.6 < a < 300 AU the M8.0-L0.5 binary frequency is 
within the range 15 ± 7%. 

Our M8.0-L0.5 binary fraction range of 15 ± 7% is 
marginally consistent with the 28 ± 9% measured for more 
massive M0-M4 dwarfs (Fischer & Marcy 1992) over the 
same separation/period range (2.6 < a < 300 AU) probed 
in this study. However, Fischer & Marcy (1992) found a bi- 
nary fraction of 32±9% over the whole range of a > 2.6AU. 
If we assume that there are no missing low mass wide bi- 
nary systems with a > 300 AU (this is a good assumption 
since such wide sep > 15" systems would have been easily 
detected in the 2MASS point source catalog as is illus- 
trated in Figure 12), then our binary fraction of 15 ± 7% 
would be valid for all a > 2.6 AU and would therefore be 
slightly lower than 32 ± 9% observed for M0-M4 dwarfs 
with a > 2.6 AU by Fischer & Marcy (1992). Hence it ap- 
pears VLM binaries (M tot < 0.185M Q ) are less common 
(significant at the 95% level) than M0-M4 binaries over 
the whole range a > 2.6 AU. 

5.2. The separation distribution function for M8.0-L0.5 

binaries 

The M8.0-L0.5 binaries are much tighter than M0-M4 
dwarfs in the distribution of their semi-major axes. The 
M8.0-L0.5 binaries appear to peak at separations ~ 4 AU 
which is significantly tighter than the broad ~ 30 AU peak 
of both the G and M star binary distributions (Duquennoy 
& Mayor 1991; Fischer & Marcy 1992). This cannot be a 
selection effect since we are highly sensitive to all M8.0- 
L0.5 binaries with sep > 20 — 300 AU (even those with 
AH > 10 mag). Therefore, we conclude that M8.0-L0.5 
stars likely have slightly lower binary fractions than G and 
early M dwarfs, but have significantly smaller semi-major 
axes on average. 

6. THE VLM BINARY POPULATION IN GENERAL 

More observations of such systems will be required to 
see if these trends for M8.0-L0.5 binaries hold over bigger 
samples. It is interesting to note that in Reid et al. (2001a) 
an HST/WFPC2 survey of 20 L stars found 4 binaries and 
a similar binary frequency of 10-20%. The widest L dwarf 
binary in Koerner et al. (1999) had a separation of only 9.2 
AU. A smaller HST survey of 10 T dwarfs by Burgasser 
et al. (2003) found two T binaries and a similar binary 
frequency of 9j^ 5 % with no systems wider than 5.2 AU. 
Therefore, it appears all M8.0-L0.5, L, and T binaries may 
have similar binary frequencies ~ 9 — 15% (for a > 3 AU). 

In Table 4 we list all the currently known VLM bina- 
ries (defined in this paper as M to t < 0.185M Q ) from the 
high-resolution studies of Basri & Martin (1999); Martin, 
Brandncr, & Basri (1999); Koerner et al. (1999); Reid et al. 
(2001a); Lane et al. (2001); Potter et al. (2002a); Burgasser 
et al. (2003); Bouy et al. (2003). As can be seen from Fig- 
ure 13 VLM binaries have a ~ 4 AU peak in their separa- 
tion distribution function with no systems wider than 15 



AU. 

From Figure 14 we see that most VLM binaries have 
nearly equal mass companions, and no system has q < 0.7. 
This VLM q distribution is different from the nearly flat q 
distribution of M0-M4 stars (Fischer & Marcy 1992). Since 
the HST surveys of Martin, Brandner, & Basri (1999); 
Reid et al. (2001a); Burgasser et al. (2003); Bouy et al. 
(2003) and our AO surveys were sensitive to 1.0 > q > 0.5 
for systems with a > 4 AU, the dearth of systems with 
0.8 > q > 0.5 in Figure 14 is likely a real characteristic 
of VLM binaries and not just a selection effect of insensi- 
tivity. However, these surveys become insensitive to tight 
(a < 4 AU) systems with q < 0.5, hence the lack of detec- 
tion of such systems may be purely due to insensitivity. 

6.1. Why are there no wide VLM binaries? 

It is curious that we were able to detect 8 systems in 
the range 0.1 — 0.25" but no systems were detected past 
0.5" (~ 16 AU). This is surprising since we (as well as the 
HST surveys of Martin, Brandner, & Basri (1999); Reid 
et al. (2001a); Burgasser et al. (2003); Bouy et al. (2003)) 
are very sensitive to any binary system with separations 

> 0.5" and yet none were found. One may worry that 
this is just a selection effect in our target list from the 
spectroscopic surveys of Gizis et al. (2000) and Cruz et 
al. (2003), since they only selected objects in the 2MASS 
point source catalog. There is a possibility that such a 
catalog would select against 0.5" — 2.0" binaries if they 
appeared extended in the 2MASS images. However, we 
found that marginally extended PSFs due to unresolved 
binaries (separation < 2") were not being classified as ex- 
tended and therefore were not removed from the 2MASS 
point source catalog. For example, Figure 12 illustrates 
that no known T-Tauri binary from list of White & Ghcz 
(2001) was removed from the 2MASS point source cat- 
alog. Although, due to the relatively poor resolution of 
2MASS (FWHM - 2 - 3"), only systems with separations 

> 3" were classified as binaries by 2MASS, all the other T- 
Tauri binaries were unresolved and mis-classified as single 
stars. In any case, we are satisfied that no "wide" (0.5" < 
separation < 2") VLM candidate systems were tagged as 
extended and removed from the 2MASS point source cata- 
log. Therefore, the lack of a detection of any system wider 
than 0.5" is not a selection effect of the initial use of the 
2MASS point source catalog for targets. 

Our observed dearth of wide systems is supported by 
the results of the HST surveys where out of 16 L and 
two T binaries, no system with a separation > 13 AU 
was detected. We find the widest M8.0-L0.5 binary is 16 
AU while the widest L dwarf binary is 13 AU (Bouy et 
al. 2003), and the widest T dwarf binary is 5.2 AU (Bur- 
gasser et al. 2003). However, M dwarf binaries just slightly 
more massive (> O.2M ) in the field (Reid & Gizis 1997a) 
and in the Hyades (Reid & Gizis 1997b) have much larger 
separations from 50-200 AU. 

In Figure 15 we plot the sum of primary and secondary 
component masses as a function of the binary separation 
for all currently known VLM binaries listed in Table 4. 
It appears that all VLM and brown dwarf binaries (open 
symbols) are much tighter than the slightly more massive 
M0-M4 binaries (solid symbols). 

If we examine more massive (SpT=A0-M5) wide bi- 
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nary systems from Close et al. (1990), we see such wide 
binaries have maximum separations best fit by a max ~ 
lOOO(M tot /O.185M ) AU in Figure 15. Any system with 
such a separation (a = a max ) will be disrupted by a differ- 
ential velocity impulse (or "escape kick") to one compo- 
nent of V esc > 0.57 km/s (see the solid upper line in Figure 
15). However, if we try to fit the VLM/brown dwarf bi- 
naries (defined here as systems with M tot < 0.185M Q ) 
we find the maximum separation is better predicted by 
a m ax VL M ~ 23.2(M tot /O.185M ) AU (lower dashed line 
in Figure 15). Since these are much smaller separa- 
tions we find that any system with such a separation 
(a = a maxvLM ) will require a larger escape velocity kick 
of V esc > 3.8 km/s to become unbound. 

To try and glean if VLM/brown dwarf binaries are re- 
ally more tightly bound, we compare their binding en- 
ergy to that of more massive binaries. Figure 16 illus- 
trates how the minimum binding energy of binaries with 
M to t < 0.185M Q is 16x harder than more massive field 
M-G binaries. In other words the widest binaries with 
Mtot < 0.185M Q appear to be 16 times "more bound" 
than the widest binaries with M tot > O.185M . 

This hardening could be a relic from dissipative 
star/disk interactions with the accretion disks around each 
star (McDonald & Clarke 1995) and/or it could be from 
an ejection event, dynamical decay, or fragmentation. We 
will briefly explore some of these possibilities in the next 
sections. 

6.1.1. Can ejection explain the lack of wide VLM 
binaries? 

Rcipurth & Clarke (2001) suggest that VLM binary sys- 
tems may have been ejected from their "mini-cluster" stel- 
lar nurseries in close triple or quadruple encounters with 
more massive objects early in their lives. Consequently, 
these ejected VLM objects are starved of accretion mate- 
rial and their growth is truncated (see a cartoon of this 
scenario in Fig. 17). Sterzik & Durisen (1998) estimate 
that around 5% of ejecta from pentuple systems are bi- 
naries. The typical ejection velocity estimated by Sterzik 
& Durisen (1998); Reipurth & Clarke (2001) is - 3 km/s 
for single objects (on average these ejected single objects 
have masses> O.185M ). Hence one might hypothesize 
tight VLM binaries of similar total mass may be ejected 
at similar velocities; however, we caution that more de- 
tailed simulations are required to estimate realistic ejec- 
tion velocities. In any case, only tightly bound binaries 
will survive the ejection process. This may explain the 
additional tightness of the VLM/brown dwarf binaries in 
Figure 15. More loosely bound VLM binaries were dis- 
solved ("kicked" past the V esc — 3.8 km/s line in Fig. 15) 
when they were ejected from their mini-cluster. Hence, 
only relatively hard a < 16 AU (V esc > 3.8 km/s) low 
mass binaries have survived until today. 

However, the ejection paradigm of Rcipurth & Clarke 
(2001) as simulated in detail by Bate et al. (2002) only pre- 
dict a VLM/brown dwarf binary fraction of < 5% which 
is below the 15 ± 7% observed for M8.0-L0.5 binaries and 
the 9^4 5 observed for T dwarf binaries (Burgasser et al. 
2003). Therefore, further simulations will be required to 
see if a larger (> 9%) binary frequency can be produced 
when low mass binaries are ejected from the mini-cluster. 



6.1.2. Is the dearth of wide low mass binaries due to 

Galactic dynamical evolution? 

Over the lifetime of a binary there will be many stochas- 
tic encounters, which will increase the potential energy 
(and therefore its separation) of the binary slowly over 
time. Eventually, these encounters may also disrupt the 
binary. In fact, this disruption timescale is roughly pro- 
portional to M to t /a for separations < 200 AU according 
to the detailed models of Weinberg et al. (1987). As we 
can see from the solid line in Figure 15, a line of constant 
Mtot/ a or constant V esc (where V esc = 0.57 km/s) can fit 
the widest A0-M0 binaries, but the same line cannot fit 
the widest of the much tighter VLM binaries (open sym- 
bols). This is quite puzzling since wide (> 200 AU) VLM 
binaries should get wider throughout their lifetime until 
they reach a = a max ~ lOOO(M tot /O.185M ) AU. By this 
point, they have reached the average minimum binding en- 
ergy {—Eund min ~ 3xl0 41 erg) of field wide binaries and 
are likely to be disrupted by a differential kick of only 0.6 
km/s. 

So why do we not observe wider VLM binaries 
if they should dynamically evolve to a = a max ~ 
1000(M tot /0.185M Q )AU - 1000AU over time? It is 
critical to note that only binaries that are wider than 
lOOO(M tot /M ) AU are wide enough to significantly evolve 
in the galactic disk according the detailed models of Wein- 
berg et al. (1987). Therefore, only binaries of initial sepa- 
ration (a ) greater than a a ~ 185 AU for M to t — O.185M 
will evolve to the minimum binding energy (which corre- 
sponds to separation of a max ~ 1000 AU) over a 12 Gyr 
lifetime. Any systems formed tighter than lOOO(M t(rf /M ) 
AU will not dynamically evolve to significantly wider sep- 
arations over time. 

So why don't we detect any 185-1000 AU VLM binaries? 
The answer may be very simple. If the formation process 
that forms VLM/brown dwarf binaries cannot produce bi- 
naries with a Q > 1000(M tot /M Q ) AU then there will be 
no significant evolution of a during the lifetime of the bi- 
naries. In other words the observed a distribution will be 
similar to the initial distribution (a a ). Since we currently 
observe no VLM systems with a > 100 AU, the observed 
separation distribution for VLM binaries is likely the same 
as their initial distribution. Figures 15 and 16 suggest 
that a Q for VLM stars is strongly truncated at ~ 16 AU. 
In addition, the observed values of V esc > 3.8 km/s and 
—Ebind > 50 x 10 41 erg are therefore likely relics of the 
formation of these VLM binaries and need to be explained 
by any successful model of star/brown dwarf formation. 

6.1.3. Could VLM binaries be the decay product of a 

dissolving mini-cluster? 

In the dynamical decay of these "mini-clusters" it is 
likely that the most massive components become a hard- 
ened binary as the lower mass objects are ejected by triple 
encounters (see bottom of Figure 17). Rcipurth & Clarke 
(2001) have hypothesized that VLM binaries may also be 
ejected from these mini-clusters since they have very low 
masses and are tight enough to remain bound after ejec- 
tion. However, it is also possible that VLM binaries are not 
ejected but instead persevere until they are the final most 
tightly bound binary remaining after the "mini-cluster" 
decays. 
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The decay of 3,4, and 5-body "mini-clusters" has been 
studied by Sterzik & Durisen (1998). The final residual bi- 
naries produced by such decays have characteristics similar 
to those observed for VLM binaries. In particular, Sterzik 
& Durisen (1998) predict that binaries composed of stars 
with masses from 0.1 — 0.2M Q would have separations of 
~ 10 AU compared to ~ 30 AU for more massive stars. 
Furthermore, they predict the mass ratio q to be in the 
range 0.6-1.0. Moreover, they predict a fairly sharp cut- 
off at <~ 10 AU for "wide" low mass binaries. All these 
predictions are roughly consistent with our observations 
of the VLM binary population. 

Since the final binary produced is typically biased to- 
wards the two or three most massive members of the 
mini-cluster (true of ~ 99 — 98% cluster decays; Sterzik 
& Durisen (1998)), the likelihood of the most massive 
objects in the mini-cluster both having masses less than 
0.090M Q is rather small. Indeed, the binary fraction of 
0.2 < M tot < 0.4M Q binaries predicted by Sterzik & 
Durisen (1998) is very small (< 1%) if each member of the 
cluster are picked randomly from the IMF. Even smaller 
binary fractions would be predicted for even lower mass 
VLM binaries (M tot < O.185M ). Therefore, only a frac- 
tion of our VLM binaries (observed to have a binary fre- 
quency of ~ 9 — 15%) are likely the residual hardened 
binaries remaining after a mini-cluster dissolves. 

Recently, Durisen, Sterzik, & Pickett (2001) have mod- 
eled a "two-step IMF" to produce a more realistic bi- 
nary population from F to early M spectral types. As 
in Sterzik & Durisen (1998) the VLM binary character- 
istics predicted in general are similar to those observed. 
In particular, this is true for their models which include 
brown dwarfs in the IMF. However, they note that mod- 
els which include brown dwarfs overestimate the num- 
ber of G-dwarf binaries (BF(1.0M Q ) ~ 75%). Moreover, 
while these models do predict some binary VLM systems 
(BF(0.1Mq < 5%) they still underestimate the frequency 
of such systems compared to the 15 ± 7% observed here 
for binaries in this mass range. 

6.1.4. Can we explain the truncation at ~ 16 AU simply 
by an initial semi-major axis distribution produced 
by fragmentation? 

One could avoid the problem of the unlikelihood of VLM 
binaries surviving the dynamics of an ejection or cluster 
decay if one simply produces VLM binaries in the same 
fashion that more massive binaries are thought to form. 
This is commonly thought to be through the fragmenta- 
tion of molecular cloud cores. A fragmenting very low 
mass molecular cloud core could directly produce a VLM 
binary without evoking any ejection or decay processes. 

It would be interesting to see if we can approximate 
a VLM binary semi-major axis distribution produced by 
fragmentation by scaling distributions of more massive bi- 
naries. If we assume M ioi /a D is roughly constant for wide 
binaries (see Figure 15), then we expect the initial a dis- 
tribution to be somewhat tighter for VLM binaries com- 
pared to more massive binaries. We can estimate the a Q 
for solar mass binaries from the well known young T-Tauri 
star a distributions. Young T-Tauri binaries are strongly 
suspected to have formed by fragmentation (White & Ghez 
2001). We can utilize a HST sample of 29 T-Tauri binary 



systems in Taurus to produce an initial separation distri- 
bution. Such a distribution from the sample of White & 
Ghez (2001) has an a peak at ~ 30 AU (as also found by 
Lcinert et al. (1993); Ghez et al. (1993)). We can scale 
this to match our observed VLM binary peak of ~ 4 AU by 
simply scaling the distribution by a OVLM <~ 0.13(a OTTau ). 
Similarly we can scale the masses of the T Tauri binaries 
to also match our mean VLM binary mass of ~ 0.15M Q 
by multiplying by the same factor of 0.13. In this manner 
we have created a plausible fragmentation-produced VLM 
a Q distribution. However, we have assumed that the other 
properties of the distribution (such as the FWHM and tail 
distributions) can be scaled along with the mean a a and 
M to t- Assuming such a scaling is valid globally, we see 
that this fragmentation a OVLM distribution has ~ 26% of 
systems wider than ~ 40 AU. Hence if our observed VLM 
distribution is to match this scaled fragmentation a Q dis- 
tribution, we would expect ~ 9 of the 34 VLM/brown 
dwarf systems observed (see Table 4) to have separations 
greater than 40 AU. However, from Figure 15 no VLM 
systems are observed to have separations greater than 16 
AU. Therefore, it appears very hard to explain the total 
lack of systems with separations greater than 16 AU by 
scaling the observed a distribution of T-Tauri stars. In 
other words, if we scale the distribution to match the VLM 
~ 4 AU peak we over-produce wide systems (26% of sys- 
tems wider than 40 AU) compared to observations. It may 
be that our "toy model" of simply scaling the T Tauri a a 
distribution is too simple an approach, but, it is a first step 
which shows that such a sharp cut-off (a a < 16 AU) is not 
easily produced by a fragmentation distribution. Detailed 
fragmentation simulations at the lowest masses will be re- 
quired to prove whether fragmentation can produce the 
VLM binary separation distribution observed. 

6.1.5. What effect does an ejection "kick" have on the 
binary distribution? 

Assuming (as has been suggested by Reid et al. (2002); 
Bate et al. (2002)) that eventually a close multiple en- 
counter ejects VLM (< 0.185M Q ) binaries, it would be 
interesting to see if such an "ejected" VLM binary dis- 
tribution would have a maximum separation of ~ 16 AU 
as observed. The full dynamical modeling of the effect 
of a binary's ejection from a mini-cluster is beyond the 
scope of this paper. However, we can simply explore, as a 
"toy-model" , the possibility that the ejection process pro- 
vides a hard differential velocity kick of order ~ 3km/s to 
the two components of an ejected VLM binary. We take 
the value of ~ 3 km/s since this is the value calculated 
for the ejection velocity of on average more massive single 
objects from a dissolving mini-cluster (Sterzik & Durisen 
1998; Reipurth & Clarke 2001). We assume, this may be 
of the right order of magnitude for the differential velocity 
kick between the VLM members of the ejected binary. We 
caution that this is only a simple toy-model, the true dif- 
ferential velocity applied during the ejection can only be 
calculated by detailed simulations of each ejection event. 

In our simple toy-model ejection only binaries with 
Vesc ^ 3.0 km/s will survive the ejection process. More- 
over, systems that are more tightly bound than V esc = 3.0 
will hardly be effected by such a kick in the tidal limit of 
Weinberg et al. (1987). To survive in the general galactic 
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field we see (from Figure 15) V esc > 0.6 km/s. Hence one 
may expect today's population of VLM binaries to all have 
VeacvLM Z (3-0 + 0.6) km/s. 

It is very interesting to note that we "observe" (from 
figure 15) V esc > 3.8 km/s for wide VLM binaries which 
is similar to the "predicted" value of V escvLM > 3.6 km/s. 
Therefore, the observed cut-off at 16 AU could be pro- 
duced by a population of VLM binaries where each system 
has been subjected to a differential kick of ~ 3 km/s in the 
ejection process. We conclude that it is possible to produce 
the observed distribution of semi-major axes for VLM bi- 
naries by applying a strong differential velocity kick to each 
binary in the distribution. However, it is still not clear if 
a VLM binary frequency of 9 — 15% could be produced by 
ejection, since ejection of a VLM tight binary may be a 
very rare event. 

7. FUTURE OBSERVATIONS 

Future observations of all of these binaries should deter- 
mine if there is still Li present in their spectrum. The most 
useful Li observation would spatially resolve both compo- 
nents in the visible spectrum. Currently there is no visible 
wavelength AO systems capable of guiding on a V ~ 20 
source and obtaining resolutions better than 0.1" in the op- 
tical (Close 2000, 2002). Therefore, we will have to carry 
out such observations from space with HST/STIS perhaps. 
Trigonometric parallax measurements should also be ob- 
tained in the near future. Within the next few years one 
should also be able to measure the masses of both compo- 
nents of several of these systems as they complete a signif- 
icant fraction of their orbit. Hence, further observations 
of these systems are a high priority for the calibration of 
the brown dwarf mass-age-luminosity relation. Continued 
searches for new VLM binaries will help define this very 
interesting and important binary population with greater 
precision. 

8. CONCLUSIONS ABOUT THE VLM POPULATION IN 
GENERAL 

Based on all 34 VLM (M tot < 0.185M©) systems cur- 
rently known from this work and that of Basri & Martin 
(1999); Martin, Brandner, & Basri (1999); Koerner et al. 
(1999); Reid et al. (2001a); Lane et al. (2001); Potter et 
al. (2002a); Burgasser et al. (2003); Bouy et al. (2003), we 
find the following general characteristics: 

• VLM binaries are tight (peak separation <~ 4 AU) 
with no systems wider than 16 AU (this is <~ 10 
times tighter than the slightly more massive M0-M4 
binaries); 

• they tend to have nearly equal mass companions 
(q <~ 0.9) with no detected companions below ~ 70% 
of the primary's mass; 

• they have a corrected binary frequency of 15±7% for 
spectral types M8-L0.5, and 9±\ 5 % for T dwarf bi- 
naries (Burgasser et al. 2003) for separations a > 2.6 
AU. This is less than the 32 ± 9% binary frequency 
of M0-M4 binaries (Fischer & Marcy 1992) and the 
~ 50% binarity of G dwarf binaries (Duquennoy & 
Mayor 1991) for similar separations. 



We find the formation of such a VLM distribution to be 
problematic with current simulations of binary star forma- 
tion. 

• Ejecting VLM binaries from their formation mini- 
clusters (Reipurth & Clarke 2001) can likely pro- 
duce the observed cut-off in separation at 16 AU, 
but ejection simulations (Bate et al. 2002) produce a 
binarity of < 5% which is smaller than the ~ 9 — 15% 
binarity observed. 

• If we ignore ejection and assume that these VLM 
binaries are instead the residual (most massive) sys- 
tem remaining from the dynamical decay of the 
mini-cluster one can likely also reproduce the tight 
cut-off at 16 AU. But the strong bias to higher 
masses in the final binary predicts that too few VLM 
binaries would be made in this fashion -assuming 
each component is picked randomly from the IMF 
(Sterzik & Durisen 1998; Durisen, Sterzik, & Pickett 
2001). 

• We briefly look at a "toy-model" of a VLM binary 
population produced by fragmentation. It appears 
that simply scaling an observed fragmentation pro- 
duced T-Tauri distribution to force a peak at 4 
AU produces systems with separations much greater 
than the 16 AU cut-off observed. 

Hence, it is not obvious that the characteristics of the 
VLM binary distribution are accurately predicted by any 
of these methods. Future detailed modeling (taking into 
account circumstellar gas disks and their effects on dy- 
namical interactions/ejections) will be required to see if 
ejection or dynamical decay can be common enough to ex- 
plain the <~ 9 — 15% binarity of VLM systems. Future 
detailed fragmentation simulations of the smallest cores 
will be needed to discern if some type of truncation occurs 
at low masses that limits VLM binaries to separations of 
less than 16 AU. 
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Table 1 

M8.0-L0.5 Stars Observed with No Likely Physical Companions Between 0.1 



2MASS 


other name 


K 3 


SpT 


Ref. 


2MASSW J0027559+221932 


LP 349-25 


9 


.56 


M8.0 


1 


2MASSW J0140026+270150 




11 


.44 


M8.5 


1 


2MASSI J0149089+295613 




11 


.99 


M9.5 


1 


2MASSW J0253202+271333 




11 


.45 


M8.0 


1 


2MASSW J0320597+ 185423 


LP412-31 


10 


.57 


M9.0 


1 


2MASSI J0335020+234235 




11 


.26 


M8.5 


1 


2MASSW J0350573+ 181806 


LP 413-53 


11 


.76 


M9.0 


1 


2MASSW J0354013+231633 




11 


.07 


M8.5 


1 


2MASSW J0810586+142039 




11 


.61 


M9.0 


1 


2MASSW J0853361-032931 


LHS 2065 


9 


.98 


M9.0 


2 


2MASSW J0928256+423054 




11 


.07 


M8.0 


3 


2MASSW J1019568+732408 




11 


.78 


M8.0 


3 


2MASSW J1124048+380805 




11 


.57 


M8.0 


3 


2MASSW J1224522-123835 


BPJ 1222-1221 


11 


.37 


M9.0 


2 


2MASSP J1309219-233035 




10 


.67 


M8.0 


4 


2MASSW J1403223+300755 




11 


.63 


M8.5 


1 


2MASSW J1421314+182740 




11 


.93 


M9.5 


1 


2MASSW J1444171+300214 


LP326-21 


10 


.57 


M8.0 


1 


2MASSW J1457396+451716 




11 


.92 


M9.0 


1 


2MASSI J1501081+225001 


TVLM 513-46546 


10 


.72 


M8.5 


5 


2MASSW J1551066+645704 




11 


,73 


M8.5 


1 


2MASSW J1553199+140033 




11 


.85 


M9.0 


1 


2MASSW J1627279+810507 




11 


.87 


M9.0 


1 


2MASSW J1635192+422305 




11 


.80 


M8.0 


1 


2MASSW J1707183+643933 




11 


.39 


M9.0 


1 


2MASSW J1733189+463359 




11 


.86 


M9.5 


1 


2MASSI J2234138+235956 




11 


.81 


M9.5 


1 


2MASSI J2334394+193304 




11 


.64 


M8.0 


1 


2MASSW J2347368+270206 




12 


.00 


M9.0 


1 


2MASSW J2349489+122438 


LP 523-55 


11 


.56 


M8.0 


1 



References. — (l)Gizis et al. (2000); (2)Kirkpatrick et al. (1995); (3)Cruz et 
al. (2003); (4)Kirkpatrick et al. (1997); (5)Tinncy et al. (1993) 
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Table 2 
The binary systems observed 



System A .7 AH AK' Sep. (") PA Age (Gyr) 



2MASSW J0746425+200032 a 





.60 ± 0.20 


0.48 ± 0.15 


0.44 ± 0. 


15 





121 ± 0.008° 


85.7 ± 1.45°° 


r n +2.5 

°- U -4.15 


2MASSW J1047127+402644 b 





.85 ± 0.25 


0.91 ± 0.20 


0.50 ± 0. 


15 





122 ± 0.008 d 


328.36 ± 3.75 od 


5.ot;-; 


LHS 2397a g 


3 


.83 ± 0.60 


3.15 ± 0.30 


2.77 ± 0. 


10 





.207 ± 0.007° 


151.98 ± 1.20°° 


7.2±« 


2MASSW J1127534+741107 





.33 ± 0.11 


0.27 ± 0.10 


0.25 ± 0. 


07 


0. 


.246 ± 0.008° 


80.23 ± 1.72°° 




2MASSW ,11311391+803222 





.13 ± 0.10 


0.15 ± 0.09 


0.14 ± 0. 


05 


0. 


.267 ± 0.006 d 


168.15 ± 0.48 od 


5.0tJ-s 


2MASSW J1426316+155701 





.78 ± 0.05 


0.70 ± 0.05 


0.65 ± 0. 


10 


0. 


.152 ± 0.006° 


344.1 ± 0.7°° 




2MASSW J2140293+162518 





.77 ± 0.05 


0.73 ± 0.04 


0.75 ± 0. 


04 


0. 


.155 ± 0.005 f 


134.30 ± 0.5 ot 




2MASSW J2206228-204705 





.17 ± 0.04 


0.08 ± 0.03 


0.08 ± 0. 


03 


0. 


168 ± 0.007 f 


68.2 ± 0.5 ot 




2MASSW J2331016-040618 


2 


78 ± 0.04 


2.64 ± 0.05 


2.44 ± 0. 


03 





573 ± 0.008 f 


302.6 ± 0.4° f 





a discovcry paper Rcid ct al. (2001a) 
b discovery paper Rcid ct al. (2002) 
Observations made on Feb 7, 2002 UT 
Observations made on April 25, 2002 UT 
"^observations made on June 20, 2001 UT 
Observations made on Sept 20, 2001 UT 
s discovery paper Freed, Close, Sz Sicglcr (2003) 
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Table 3 

Summary of the new binaries' A & B components 



Name J H Ks M Ks c SpT a Est. Mass b Est. D (pc)° Sep. (AU) P (yr) d 



2M0746A 


12 


.23 


± 


0. 


083 


11. 


.53 


± 





07 


11. 


04 


± 





07 


10.58 ± 0.07 


L0.5 


0. 


082+.° 


001 

.004 


12 


.34 ± 0.05° 


1. 


.47 ± 0.11 


11 


+3 

-2 


2M0746B 


12 


.83 


± 


0. 


21 


12 


.01 


± 





16 


11. 


48 


± 





16 


11.0 ± 0.17 


L2 


0. 


.078+.° 


001 
009 














2M1047A 


12 


.85 


± 


0. 


08 


12 


.00 


± 


0. 


.07 


11. 


80 


± 





.07 


9.95 ± 0.39 


M8 





.092+° 


009 
012 


23 


.4 ± 4.3 


2. 


.9 ± 0.6 


11 


+ 5 
-3 


2M1047B 


13 


.70 


± 


0. 


26 


13 


.00 


± 


0. 


21 


12. 


31 


± 





.16 


10.46 ± 0.42 


L0 


0. 


084+.°, 


(106 
016 














LHS 2397aA 


11 


.86 


± 


0. 


05 


11 


.32 


± 


0. 


.05 


10. 


.80 


± 





.03 


10.03 ± 0.09 


M8 





■090±° 


(104 


14 


.3 ± 0.4° 


3. 


86 ± 0.18 


22 


+3 
-3 


LHS 2397aB 


15 


69 


± 


0. 


60 


14 


.47 


± 


0. 


30 


13. 


.57 


± 





.10 


12.80 ± 0.12 


L7.5 





.068jlg 


001 
007 














2M1127A 


13 


.66 


± 


0. 


.06 


12 


.00 


± 


0. 


.05 


12. 


.60 


± 





.05 


9.96 ± 0.37 


M8 





■092±° 


009 
.012 


33 


.67 ±5.9 


8. 


31 ± 1.49 


57 


+ 26 
-15 


2M1127B 


13 


.99 


± 


0. 


.13 


13 


.26 


± 


0. 


11 


12. 


.85 


± 





.08 


10.21 ± 0.38 


M9 





.087+.° 


007 
.013 














2M1311A 


13 


.40 


± 


0. 


06 


12 


.82 


± 


0. 


.05 


12. 


39 


± 





.04 


10.07 ± 0.37 


M8.5 





.089+.° 


008 
.011 


29 


.1 ± 5.1 


7. 


7± 1.3 


51 


+ 24 
-13 


2M1311B 


13 


.62 


± 


0. 


.11 


12 


.97 


± 


0. 


10 


12. 


53 


± 





.07 


10.21 ± 0.38 


M9 





.087±° 


007 
.013 














2M1426A 


13 


.30 


± 


0. 


.04 


12 


.63 


± 


0. 


04 


12. 


18 


± 





.05 


10.09 ± 0.37 


M8.5 





.088+0 


009 
.010 


26 


.1 ± 4.5 


3 


.9 ± 0.7 


10 


+9 
-5 


2M1426B 


14 


.08 


± 





.06 


13 


.33 


± 


0. 


06 


12. 


.83 


± 





.11 


10.74 ± 0.38 


LI 





o7 6 ;° 


008 
.014 














2M2140A 


13 


.37 


± 





.04 


12 


.71 


± 


0. 


.04 


12. 


22 


± 





.04 


10.17 ± 0.37 


M9 





.088+.° 


007 
.012 


25 


.6 ± 4.4 


3 


.9 ± 0.6 


10 


+ 9 
— 5 


2M2140B 


14 


14 


± 





06 


13 


.44 


± 


0. 


.05 


12. 


.97 


± 





.05 


10.92 ± 0.37 


L2 





.078+.°, 


004 
020 














2M2206A 


13 


.10 


± 





.04 


12 


.46 


± 





04 


12. 


.06 


± 





.04 


9.93 ± 0.37 


M8 





.092 + o 


.000 
.011 


26 


.7 ± 4.5 


4. 


4 ± 0.7 


22 


+ 10 

— 5 


2M2206B 


13 


.27 


± 





.05 


12 


.54 


± 





05 


12. 


.14 


± 





.05 


10.01 ± 0.37 


M8 





■091±° 


008 
.011 














2M2331A 


13 


.02 


± 





.04 


12 


.38 


± 





04 


12. 


.03 


± 





.04 


9.93 ± 0.08 


M8.0 





■093+.° 


002 
.002 


26 


.2 ± 0.6° 


15.0 ± 0.37 


159 -39 


2M2331B 


15 


.80 


± 


0. 


05 


15 


.02 


± 


0. 


06 


14. 


.47 


± 





.05 


12.37 ± 0.12 


L7 





.067+.° 


002 
013 















a Spcctral types estimated by 3.97xM/^ s — 31.67 (Dahn et al. 2002; Sicglcr et al. 2003) with ±1.5 spectral subclasses of error in these estimates 
(note a value of SpT— 10 is defined as L0 in the above equation). 

b Masses (in solar units) from the models of Chabrier et al. (2000) —sec Figure 3. 

c Photometric distances estimated for the primaries by M Ks — 7.71 + 2.14( J — Ks) which is valid for Ml < SpT < LI (Siegler et al. 2003). 

d Pcriods estimated assuming face-on circular orbits. In the cases of 2M0746 and LHS 2397a, the larger of the 2 observed separations (2.7 AU 
& 3.86 AU from Reid et al. (2001a) and Freed, Close, & Siegler (2003), respectively) have been used to more accurately estimate the period. 

e These systems have trigonometric parallaxes. 
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Table 4 

All Known Resolved VLM Binaries 



Name 


Sop. 


Est. 


Est. M A 


Est. M B 


Est. 


Period 13 


Rof.° 




AU 


SpT A /SpT B 


M Q 


M e 








PPL 15 a 


0.03 


M7/M8 


0.07 


0.06 


5.8 


days 





Gl 569B 


1.0 


M8.5/M9.0 


0.063 


0.06 


3 




1,2 


SDSS 2335583-001304 


1.1? 


L17/L4? 


0.079 


0.074 


3 




12 


2MASSW J1112256+354813 


1.5 


L4/L6 


0.073 


0.070 


5 




12 


ZiVl/Yool .J lo.y41:Jo-zuozz 1 


1 Q 


1 O.O/ 1 o.o 


v.vo 


v.vo 


Q 
O 




y 


2MASSW J0856479+223518 


2.0 


L57/L8? 


0.071 


0.064 


8 




12 


DENIS-P J185950. 9-370632 


2.0 


L0/L3 


0.084 


0.076 


7 




12 


HD130948B 


2.4 


L2/L2 


0.07 


0.06 


10 




3 


2MASSW J0746425+200032 


2.7 


L0.5/L2 


0.082 


0.078 


12 




4 & this paper 


2MASSW J1047127+402644 


2.7 


M8/L0 


0.092 


0.084 


11 




8 & this paper 


DENIS-P J035726. 9-441730 


2.8 


L2/L4 


0.078 


0.074 


12 




12 


2MASSW J0920122+351742 


3.2 


L6.5/L7 


0.068 


0.068 


16 




4 


LP415-20 


3.5 


M7/M9.5 


0.095 


0.079 


15 




11, this survey 


2MASSW J1728114+394859 


3.7 


L7/L8 


0.069 


0.066 


19 




12 


LHS 2397a 


3.9 


M8/L7.5 


0.090 


0.068 


22 




7, this paper 


2MASSW J1426316+155701 


3.9 


M8.5/L1 


0.088 


0.076 


19 




this paper 


2MASSW J2140293+162518 


3.9 


M9/L2 


0.092 


0.078 


22 




this paper 


2MASSW J2206228-204705 


4.4 


M8/M8 


0.092 


0.092 


22 




this paper 


2MASSs J0850359+105716 


4.4 


L6/L8 


0.05 


0.04 


30 




4 


2MASSW J1750129+442404 


4.8 


M7.5/L0 


0.095 


0.084 


25 




11, this survey 


DENIS-P J1228. 2-1547 


4.9 


L5/L5 


0.05 


0.05 


34 




10 


2MASSW J1600054+170832 


5.0 


L1/L3 


0.078 


0.075 


29 




12 


2MASSW J1239272+551537 


5.1 


L5/L5 


0.071 


0.071 


31 




12 


2MASSI J1553022+153236 


5.2 


T7/T7.5 


0.040 


0.035 


43 




9 


2MASSW J1146345+223053 


7.6 


L3/L4 


0.055 


0.055 


63 




6 


2MASSW J1311391+803222 


7.7 


M8.5/M9 


0.089 


0.087 


51 




this paper 


2MASSW J1127534+741107 


8.3 


M8/M9 


0.092 


0.087 


57 




this paper 


LP475-855 


8.3 


M7.5/M9.5 


0.091 


0.080 


58 




11, this survey 


DENIS-P J0205. 4-1159 


9.2 


L7/L7 


0.07 


0.07 


75 




6 


2MASSW J2101349+175611 


9.6 


L7/L8 


0.068 


0.065 


82 




12 


2MASSW J2147436+143131 


10.4 


L0/L2 


0.084 


0.078 


83 




12 


2MASSW J1449378+235537 


11.7 


L0/L3 


0.084 


0.075 


100 




12 


DENIS-P J144137. 3-094559 


13.5 


Ll/Ll 


0.079 


0.079 


124 




12 


2MASSW J2331016-040618 


15.0 


M8.0/L7 


0.093 


0.067 


159 




this paper 



a PPL 15 is a spectroscopic binary (Basri & Martin 1999) 

b This "period" is simply an estimate assuming a face-on circular orbit 

c REFERENCES-(0)Basri Sz Martin (1999); (1) Kenworthy et al. (2001); (2) Lane ct al. (2001); (3) Potter ct al. 
(2002a); (4)Reid et al. (2001a); (5) Martin, Brandncr, & Basri (1999); (6) Koerner et al. (1999); (7) Freed, Close, 
& Sicgler (2003); (8) Rcid ct al. (2002); (9) Burgasscr ct al. (2003); (10) Martm, Brandner, & Basri (1999); (11) 
Sicglcr ct al. (2003); (12) Bouy et al. (2003) 

Gl 569B and HD130948B arc binary brown dwarfs that orbit normal stars, for AO observations these bright 
primary stars were guided on, not the brown dwarfs 

G We define VLM binaries as 2 star systems where M to t < O.185M0. Very young evolving systems (like GG 
TauBaBb (White & Ghez 2001)) are not included, nor arc over-luminous systems which are not resolved into 
binaries. 
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2M0746 — \ 


2M1 047^— ~\ 










0.1 1 9" 


2M1127_^^ 


2M1311 


0.246" 





Fig. 1.— The 12x10s X' image of the L0.5 2MASSW J0746425+200032 binary (which is the one overlap object between our study and 
that of Reid et al. (2001a)). Note the very different position (PA = 85.7°, sep = 0.121" = 1.49 AU on 2002/2/7) of our image compared 
to the PA=168.8° sep 2.7 AU published by Bouy et al. (2003) from HST observations of Reid et al. (2001a) at 2000/4/15. This suggests 
that the orbit could have a period of ~ 4yr. In the near future this system could have an orbital solution. At a resolution of < 0.1" both 
components are clearly visible. We also show K' images of the new binaries 2MASSW J1047127+402644, 2MASSW J1127534+741107, & 
2MASSW J1311391+803222. The pixels are 0.199"pix~ 1 . The contours are linear at the 90, 75, 60, 45, 30, 15, and 1% levels. North is up 
and east left in each image. 
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Fig. 2. — In figure (a) we see the 12x10s K' image of the 2MASSJ 1426316+155701 binary discussed in Close et al. (2002b) at a resolution 
of 0.131". In Figures (b-d) we show K' images of the binaries 2MASSW J2140293+162518, 2MASSWJ 2206228-204705, and 2MASSWJ 
2331016-040618, respectively. The pixels are 0.0199" pi* -1 - The contours are linear at the 90, 75, 60, 45, 30, 15, and 1% levels. North is up 
and east left in each image. 
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0.04 0.00 0.08 0.10 

M./M 

Fig. 3. — The latest Chabrier et al. (2000) DUSTY evolutionary models arc shown custom integrated over the Ks bandpass. The locations 
of the 2 components of 2M0746 arc indicated by the crosses. The polygon encloses the error in the Mjf s and age range 0.6 — 7.5 Gyr of the 
system. The Mk s of each secondary is determined by the addition of AKs plus the Mk s of the primary. The errors for the primary are 
enclosed by the upper polygon (outlined by a solid thick line), the secondary is bounded by the lower polygon (outlined by a thick dashed 
line). The models suggest a primary mass of 0.082 M Q with a range 0.078 - O.O83M with temperatures of 2375 K (2323-2409K). For the 
secondary the models suggest a mass of 0.078Mq with a range 0.068 — O.O79M0 with temperatures of 2173K (2034-2252 K). The isochrones 
plotted are 0.3, 0.6, 0.65, 0.7, 0.85, 1.2, 1.7, 3.0, 5.0, 7.5, & 10.0 Gyr. The isotherms run in equal intervals from 2600 K to 1400 K in steps of 
200 K. 
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0.04 0.00 0.08 0.10 

M./M 

Fig. 4. — The same as figure 3 except for 2M1047. Note the larger errors in Mks compared to 2M0746 since we have to currently rely on a 
photometric parallax for this system. The models suggest a primary mass of 0.092 Mq with a range 0.079 — O.1O2M0 with temperatures of 
2660 K (2460-2830K). For the secondary the models suggest a mass of 0.084M Q with a range 0.068 - O.O9OM with temperatures of 2430K 
(2164-2625 K). 
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Fig. 5. — Figure 3 from Freed, Close, & Siegler (2003). Note that the isochrones are different for this plot compared to all the other plots 
in this paper (but the isotherms are the same). More detail about LHS 2397a can be found in Freed, Close, & Siegler (2003). LHS 2397aB 
(mass 0.068 Mq) is one of the tightest brown dwarf companions ever to be imaged around a star. 
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Fig. 6. — The same as for figure 3 except for 2M1127. Note the larger errors in Mks compared to 2M0746 since we have to currently rely 
on a photometric parallax for this system. The models suggest a primary mass of 0.092 Mq with a range 0.080 — 0.101M© with temperatures 
of 2650 K (2460-2810K). For the secondary the models suggest a mass of O.O87M with a range 0.074 - O.O95M with temperatures of 2540K 
(2330-2710 K). 
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Fig. 7. — The same as for figure 3 except for 2M1311. Note the larger errors in Mks compared to 2M0746 since we have to currently rely 
on a photometric parallax for this system. The models suggest a primary mass of 0.089 Mq with a range 0.078 — 0.098Mq with temperatures 
of 2610 K (2400-2760K). For the secondary the models suggest a mass of O.O88M with a range 0.074 - O.O95M with temperatures of 2540K 
(2331-2710 K). 
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Fig. 8. — The same as for figure 3 except for 2M1426. Note the larger errors in Mk s compared to 2M0746 since we have to currently rely 
on a photometric parallax for this system. The models suggest a primary mass of 0.088 Mq with a range 0.077 — 0.097Mq with temperatures 
of 2560 K (2400-2750K). For the secondary the models suggest a mass of 0.076M Q with a range 0.062 - O.O85M with temperatures of 2280K 
(2016-2480 K). 
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Fig. 9. — Same as for figure 3 except for 2M2140. Note the larger errors in Mk s compared to 2M0746 since we have to currently rely on a 
photometric parallax for this system. The models suggest a primary mass of 0.088 Mq with a range 0.075 — O.O95M0 with temperatures of 
2560 K (2360-2720K). For the secondary the models suggest a mass of O.O78M with a range 0.057 - 0.083M Q with temperatures of 2210K 
(1880-2390 K). 
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Fig. 10. — The same as for figure 3 except for 2M2206. Note the larger errors in Mks compared to 2M0746 since we have to currently rely 
on a photometric parallax for this system. The models suggest a primary mass of 0.092 Mq with a range 0.081 — 0.102Mq with temperatures 
of 2670 K (2480-2830K). For the secondary the models suggest a mass of O.O91M with a range 0.079 - 0.100M Q with temperatures of 2640K 
(2440-2790 K). 
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Fig. 11. — The same as for figure 3 except for 2M2331. This system is associated with the F8 star HD221356 (Gizis et al. 2000) which has a 
Hipparcos parallax of 26.2 ± 0.6 pc. The system is very wide with a separation of 0.057 pc between HD221356 and the 2M2331 system (Gizis 
et al. 2000). The models suggest a mass for 2M2331A of 0.093 M with a range 0.091 - 0.095M Q with temperatures of 2670 K (2640-2700K). 
For 2M2331B the models suggest a mass of 0.067M Q with a range 0.055 - O.O7OM with temperatures of 1560K (1470-1584 K), making 
2M2331B a definite brown dwarf companion. 
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Fig. 12. — It is possible that the reason none of these 2MASS based surveys have detected wide low mass binaries is because the 2MASS 
pipeline flags blended binaries as "extended" and they are left out of the point source catalog completely. To check the validity of this 
possibility we plot above all the T Tauri binaries from White & Ghcz (2001) found in the 2MASS second release point source catalog. None of 
these binaries were classified as extended sources sources, hence they each had at least one entry in the 2MASS point source catalog. However, 
the low angular resolution of 2MASS resulted in only systems with separations greater than ~ 4" actually being classified as binaries. All 
systems tighter than 3" were unresolved and classified as single stars. A fainter sample of white dwarf binaries yielded the same result. We 
found no examples of binaries tighter than 2" being removed from the point source catalog. Therefore, we feel the lack of detection of wide 
low mass systems is not a selection effect from the initial use of the 2MASS point source catalog to define candidate low mass systems by 
Gizis et al. (2000); Kirkpatrick ct al. (2000); Burgasser et al. (2003); Cruz et al. (2003). 
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Fig. 13. — Here we plot all 34 VLM binaries from Table 4. Note how VLM binaries appear to have smaller separations compared to the 
M0-M4 binaries of Fischer & Marcy (1992). Both distributions are normalized to unity binary fraction. We have not tried to correct for 
instrumental incompleteness. Hence we underestimate the number of VLM binaries with a < 2 AU. Poisson error bars are plotted; the sharp 
peak at 4 AU and the lack of any wide (a > 16AU) systems are real features of the distribution and arc significantly different from that 
observed in more massive M0-M4 binaries (Mtot ^ 0.3Mq). Binary systems with Mtot < 0.185Mq appear ~ 10 times tighter compared to 
just slightly (2 — 5 times) more massive M0-M4 binaries. We have converted the semi-major axis distribution values of Fischer & Marcy 
(1992) to observed separation by dividing the semi-major axis values by 1.26 (sec equation 7 in Fischer & Marcy (1992)). 
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Fig. 14. — Here we plot the 20 widest (a > 3.8 AU) VLM binaries from Table 4. We see that VLM binaries appear to have a q mass ratio 
peaked near unity. This is quite different from the much flatter distribution of M0-M4 binaries (dotted line) of Fischer & Marcy (1992). The 
difference is likely real, and not a insensitivity effect, since the high-resolution AO and HST observations should be sensitive to mass ratios 
as low as q ~ 0.5 if a > 3.8 AU (Bouy et al. 2003; Burgasser et al. 2003). However, for very faint companions both HST/WFPC2 and AO 
surveys are likely incomplete for tight (a < 8 AU) VLM binaries with q < 0.5. Therefore, no conclusions can currently be drawn about the 
likelihood of VLM systems with q < 0.5 since we are largely insensitive to such systems in the current studies. 
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Fig. 15. — Here we plot all 34 currently known VLM binaries. The nine M8.0-L0.5 systems of this paper plus seven more VLM M binaries 
from Table 4 are plotted as open six sided stars. The L dwarf binaries of Koerner et al. (1999); Martin, Brandner, & Basri (1999); Reid 
et al. (2001a); Bouy et al. (2003) are shown as open triangles, and the two T dwarf binaries of Burgasser ct al. (2003) as open four-sided 
stars. For comparison, we have plotted all the visual M0-A0 binaries inside 25pc from Close et al. (1990) as solid triangles. In addition, low 
mass Hyades binaries from Reid & Gizis (1997b) are plotted as solid circles, low mass field M-dwarfs from Rcid & Gizis (1997a) are plotted 
as solid pentagons. All A0-M5 star/brown dwarf systems arc plotted as solid squares (Reid et al. 2001a). Note how there are no low mass 
systems with separations > 16 AU. It is seen that for more massive binaries (Mtot > 0.185Mq) the maximum observed separation can be 
fit to a ma x = lOOO(Mt o t/O.185M0) AU (solid line). When a = a max the escape velocity is V eS c = 0.57 km/s. However, for VLM binaries 
(where Mtot < 0.185MQ) the best fit of the widest systems is a maXvLM = 23.2(Mt o t/0.185MQ) AU (dashed line). For VLM systems with 
a = a ma x VLM the escape velocity is 3.8 km/s. Hence VLM binaries appear sharply tighter (a < 16 AU) and have escape velocities at least 3 
km/s higher than more massive wide binaries. 
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Fig. 16. — To see if low mass binaries really are more bound (i.e. "harder") wc plot here the binding energy for all the systems in Figure 
15. The data (individual masses and separations) are from Table 4 and the symbols are the same as in Figure 15. We see that the widest 
VLM binaries (open symbols) are ~ 16 times harder than their more massive wide binary counterparts (solid symbols). This effect was also 
noted by Burgasser et al. (2002). 
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Fig. 17. — A simple cartoon of the hypothesized ejection of the low mass members of a formation mini-cluster. Typically a dissolving 
mini-cluster will evolve to a final hardened binary containing the most massive stars in the cluster. The other lower mass members are ejected 
through close encounters with the more massive members. If it is possible to eject VLM binaries, only tightly bound VLM binaries could 
survive the ejection process (Rcipurth & Clarke 2001). 



